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ABSTRAGT 

We present the results of our extensive search using the Bayesian block method for long tails 
following short bursts from a magnetar, SGR J1550—5418, over all RXTE observations of the source. 
We identified four bursts with extended tails, most of which occurred during its 2009 burst active 
episode. The durations of tails range between ^13 s and over 3 ks, which are much longer than the 
typical duration of bursts. We performed detailed spectral and temporal analysis of the burst tails. 
We find that the spectra of three tails show a thermal nature with a trend of cooling throughout the 
tail. We compare the results of our investigations with the properties of four other extended tails 
detected from SGR 1900+14 and SGR 1806—20 and suggest a scenario for the origin of the tail in the 
framework of the magnetar model. 

Subject headings: pulsars: individual (SGR J1550—5418, IE 1547.0-5408, PSR J1550-5418) — stars: 
neutron — X-rays: stars 


1. INTRODUCTION 

Magnetars — neutron stars powered by their ex¬ 
tremely strong magnetic fields (|Duncan &: Thompson! 
ll992l : [^omDSon fc Duncanll 199511 — are distinguished by 
the emission of energetic bursts observed in the hard 
X-ray/soft gamma-ray band. Gurrently there are 28 
sources classified as magnetars (see the magnetar cata¬ 
log], (|01ausen &: Kaspl 120141) for detailed information). 
Magnetar bursts can be classified according to durations 
and energetics: short bursts last a fraction of a second 
and involve an isotropic energy of < 10^*^ erg. Interme¬ 
diate events are slightly longer, typically a few seconds, 
and the emitted energy is about 2 orders of magnitude 
larger. Magnetars emit giant flares but very rarely; only 
three such flares have been observed to date. The gi¬ 
ant flares are at the extreme of the burst energy scale 
(>10"^'^ erg) and relatively long, lasting a few hundreds of 
seconds, during which there are remarkable spectral and 
temporal variations. For a comprehensive list of studies 
on magnetar bursts, see the magnetar burst librarjQ. 

SGR J1550-5418, also known as IE 1547.0-5408 or 
PSR J1550-5418, is a mag netar with current ly the short¬ 
est spin period, 2.072 s (jEnoto et al.l l201(ill . The spin 
period and s pin-down rate were measured first in the 
radio band (jGamilo et al.l 120071 1. It was first pro¬ 
posed as a magnetar candidate by iGelfand &: GaensI^ 
(|2007ll based on its magnetar-like X-ray spectrum and 
association with a supernova remnant. Identifica¬ 
tion of its spin period and spin-down rate, which 
implies a magnetic field strength of 2.2x10^^ G, 
furt her supported the suggested magnetar hypothe¬ 
sis (iGamilo et al.l 1200711 . Although there were im¬ 
plications tha t it has gone through X-ray brighten¬ 
ing e pisodes (|Gelfand fc Gaenslerl l2007t iHaloern et al.l 
1200811 , magnetar-like bursts f rom the source were not ob¬ 
served until 2008 October (jKrimm et al.l [200^ . SGR 
J1550—5418 exhibited other intense bursting episodes in 

sinemsmus@sabanciuniv.edu 

^ http://www.physics.mcgill.ca/-pulsar/magnetar/main.html 

^ https://staff.fnwi.uva.nl/a.1.watts/magnetar/mb.html 


2009 January and March (iConnaughton fc Briggsl 120091 

Ivon Kienlin fc Gonnaughtonll2009lL 

There have been numerous extensive investigations in 
order to understand the burst and pe rsistent X-ray emis - 
sion properties of SGR J1550—5418. Ilsrael et all (j2010ll . 
using Swift observations of the 2008 October burst ac¬ 
tivation, found that the 2—10 keV flux was elevated by 
^50 times above its quiescent level, and that its pulsed 
fraction has also increased significantly. Spectral analy¬ 
sis of the bursts observed with the Burst Alert Telescope 
(BAT) on board Swift revealed that a blackbody model 
with temp erature of 11 keV re presents the burst s pectra 
very well ([Israel et al.ll2010ll . | von Kienlin et akl (|2012[ 1 
analyzed the bursts observed with the Fermi Gamma- 
ray Burst Monitor (GBM) in 2008 October and 2009 
March-April. They reported that the spectral character¬ 
istics of the bursts observed in two active episodes sepa¬ 
rated by about 5 months are different: the 2008 October 
burst spectra are best described with a single-blackbody 
function, while bursts observed in 2009 March-April are 
better fit with an optically thin thermal bremsstrahlung 
model. They interpreted this variation as a reflection 
of the changes in magnetic field structure of the source 
due possibly to another extreme-intense bursting episode 
that o ccurred in between these tw o periods (in January 
2009). Ivan der Horst et al.l (l2012fl analyzed Fermi GBM 
observations of 286 bursts detected during a week fol¬ 
lowing 2009 January 22, the most burst-active episode 
of the source. They reported that burst spectra can 
be described equally well w ith a Comptoni zed model 
or double-blackbody model. iLin et akl (l2012tl addition¬ 
ally used simultaneous Swift observations to analyze the 
bursts of the same period including the soft X-ray band 
and found that the double-blackbody model represents 
the spectra better than the Comptonized model. 

Besides typical short bursts, there have been reports 
on more energetic events from SGR J1 550—5418, mostly 
during its 2009 January active phase. iMereghetti et al.1 
(|2009ll reported that some energetic bursts (with ener¬ 
gies as high as 10^^ erg) detected with INTEGRAL are 
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followed by long emission episodes which are modulated 
with t he spin period of the neutron star. iKaneko et al.l 
(|2010r) identified a 150 s long enhanced persistent emis¬ 
sion phase du ring which pulsed sig nals were detected up 
to '^110 keV. iKiiiner et al.l (j2012h identified two events 
in the Rossi X-ray Timing Explorer (RXTE) observa¬ 
tions of the same active episode and called them ‘mini 
outbursts’ due to their long emission (^450 s and ~130 
s) at a much lower intensity than the bursts, but clearly 
above the persistent emission level. Here, we call the 
latter events bursts with extended tails. 

SGR 1900-1-14 and SGR 1806—20 have shown energetic 
bursts with extended tails; extended tails typically last a 
few hundreds of seconds but can be as high as thousands 
of seconds. In all of these cases, the spectral proper¬ 
ties of the tail emission are different than from those of 
the bursts, as well as those of the persistent emission; 
the tail spectra are well fitted with a blackbody model 
with decreasing temperature throughout the course of 
the tail, which implies a cooling thermal component 
on the surface , possibly heated by the in i tiating burst 
(llbrahim et al.l [ 20011 iLenters et all 120031: IGogfis et al.l 


1201111 . ILenters et al.l ( 200311 and iGogfis et al’l ( 201111 

showed that the total energy contained in the extended 
tails accounts for a constant percentage of the initiating 
burst event; it is ^2% for SGR 1900-1-14 (ILenters et al’l 
1200311 while for the two detected ta ils in SGR 1806—2 0 
the ratios are ~0.34% and ~0.63% (|Gogus et al.ll20lill . 

Motivated by the detection of extended tails from other 
magnetars, and having alread y identified two bur sts with 
tails from SGR J1550—5418 ( Kuiper et al.ll20lll , we ex¬ 
tensively searched for extended tails following bursts in 
all available RXTE observations of SGR J1550—5418. 
However, detection of these tails on short timescales is 
not optimal due to variations of the background emis¬ 
sion and, sometimes, the existence of hundreds of bursts 
in the active episode. Thus, in the work p resented here, 
we a pplied a Bayesian block algorithm (jScargle et al.l 
1201311 . which can detect local variabilities more robustly 
to search for burst tails from SGR J1550—5418. In the 
following, the RXTE observation details are found in Sec¬ 
tion [2l and in Section |3] we explain our search methodol¬ 
ogy for the detection of tails, as well as our search results. 
We present the results of our detailed spectral and tem¬ 
poral investigations of the identified burst tails in Sec¬ 
tions [4] and O respectively. We discuss the physical im¬ 
plications of our results and compare with the properties 
of extended burst tails observed from other magnetars in 
Section 0 


2. OBSERVATIONS AND DATA ANALYSIS 

RXTE, which was operational from 1995 December to 
2012 January, had two instruments on board: Propor¬ 
tional Counter Array (PGA), detecting photons in the 
energy range 2—60 keV, and The High Energy X-ray 
Timing Experiment (HEXTE) covering an energy range 
of 15—250 keV. PGA had five proportional counter units 
(PCU) labeled from 0 to 4, each consisting of one propane 
veto, three xenon, and one xenon veto layer. In the 2—10 
keV energy range, PCU sensitivity limit was 4x 10“^^ erg 
s“^ cm~ , and telemetr y rate as high as 20,000 counts 
s“^ (|Jahoda et al.ll200^ . 

Here we used all available archival data obtained bi¬ 
weekly between 2008 October and 2010 December (191 


TABLE 1 

Details of Observations That Contain Extended Tails. 


Event 

Date 

Time (UTC)** 

Active PCUs 

A 

2009 Jan 22 

22:48:45.44 

2 

B 

2009 Feb 06 

18:29:03.15 

2, 3 

C 

2009 Mar 30 

14:13:06.20 

1, 2 

D 

2010 Jan 11 

21:12:23.40 

2 


^ Denotes the start of the event. 


observations, total exposure of ~702 ks; see Figure [T] 
for the time distribution of these observations). In our 
investigations, we used data collected with the PGA 
only. For timing analysis, we converted the arrival times 
to the time at the Solar System barycenter using the 
source coordinates of R.A. = 15^50” ’' 54^.11 and deck = 
—54°18'23".7 given bv ICamilo et al.l (|2007f) . 

3. SEARCH FOR EXTENDED TAILS 

We implemented a Bayesian blocks based algorithm to 
identify extended burst tails, which immediately follow 
the bursts with count rate much lower than that of the 
burst, but still higher than the count rate of the pre-burst 
data. The Bayesian block algorithm is a segmentation 
method in order to detect local variabilities in time series 
data by separating the data into blocks of statistically 
significant variation s that maximize the likelihood value 
(|Scargle et al.ll201^ . This method was use d to search for 
weak magnetar bursts bv ILin et al.l (I2013H and resulted 
in successful detection of the dimmest bursts observed 
from magnetars. 

To optimize the extensive search, we took a two-step 
approach: a signal-to-noise ratio (S/N) burst identifica¬ 
tion followed by a Bayesian block tail search, which we 
describe below. First, we searched for burst candidates 
in the 0.125 s binned light curves of the source using an 
S/N criterion with significance 5.5cr above background in 
the 2—20 keV energy range. The S/N search identified 
878 burst candidates in total. Then, for each of the burst 
candidates we regenerated the light curves with 1-s res¬ 
olution, for the interval between 200 s before and 1000 
s after the burst time. We then applied the Bayesian 
block algorithm to the 1 s light curves and obtained the 
Bayesian block representation of the light curves. Using 
those, we searched for the increase in the count rate that 
matches the burst times found by the S/N search. Once 
found, the beginning of the identified block was taken as 
the event start time, and the count rate of the preceding 
block was assigned as the background level. Then, the 
first decrease in the count rate after the event start time 
below the assigned background level determines the pu¬ 
tative end of the tail. If the algorithm cannot determine 
the end of the tail within the selected data segment (from 
T—200 to T-l-1000 s), we extend the post-burst interval 
by 200 s and perform the search again. For our analysis 
we run two iterations, and if the algorithm still cannot 
assign an end to the tail, we conclude that there is no 
tail associated with that particular burst. We note that 
the assigned background level can be affected by the ex¬ 
istence of data gaps, and in such cases the tail may have 
been missed in the search. Therefore, to account for such 
cases, we also artificially elevated the background level to 
30% of the burst candidate’s count rate and performed 
the search again. The tail end time in such a case was the 
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Fig. 1.— Distribution of SGR J1550—5418 observations performed by RXTE. Downward arrows indicate the tails detected with RXTE 
in events A, B, C, and D. Exposures of these observations range from 0.5 to 17 ks, with a mean of about 3.7 ks. 


TABLE 2 

Duration and count rates of bursts and tails in 2—25 keV 

BAND. 


Event 

Duration (s) 

Counts s ^ ^ 

A 



Burst 

1.83 

5573 (5524) 

Tail 

17.00 

102 (53) 

B 



Burst 

0.12 

767 (737) 

Tail 1 

47.07 

68 (28) 

Tail 2 

146.91 

47 (7) 

C 



Burst 1 

0.35 

18450 (18400) 

Burst 2 

0.17 

2206 (2147) 

Burst 3 

0.20 

1270 (1218) 

Tail 1 

148.22 

74 (19) 

Tail 2 

187.83 

63 (8) 

Tail 3 

279.80 

57 (2) 

D 



Burst 1 

0.30 

203 (190) 

Burst 2 

0.45 

340 (329) 

Tail 1 

113.30 

31 (14) 

Tail 2 

647.40 

20 (3) 


® The values In parentheses are 
background-subtracted count rates. 


end of the observation. We eliminated the tail candidates 
that lasted less than 1 s since extended tails of our in¬ 
terest have much longer durations. Finally, we excluded 
false detections that are due to known data anomalies. 

Based on these criteria, we identihed a total of four 
bursts with extended tails from SGR J1550—5418, which 
we labeled A—D in chronological order for convenience 
in this paper. In Table [U we provide observational de¬ 
tails of these events. Note that our sample includes the 
two event s (events C and D) v isually identified and men¬ 
tioned in IKuioer et al.l (|2012ll . In Figure [21 we present 
the light curves of the bursts with extended tails, as well 
as their Bayesian block representations. Among these 
events, the shortest event duration determined by the 
algorithm is 15 s, and longest duration is 3534 s. For the 
longest event, the tail end time corresponds to the end 
of the observation (see Figure [2] top right panel). 

4. SPECTRAL ANALYSIS OF TAILS AND ASSOCIATED 
BURSTS 


We investigated spectral properties of the extended 
tails and the bursts that are associated with these tails 
in order to uncover their emission properties and com¬ 
pare them with the detected tails from SGR 1900+14 
and SGR 1806-20. We used HEASOFT v6.16 to per¬ 
form data extraction. First, we applied standard filter¬ 
ing to the PGA data (Earth occultations. South Atlantic 
Anomaly passages, electron contamination, etc.). We 
perfor med spectral fi ts using XSPEG software version 
12.8.2 (jArnaudlll996[l in energy range 2—25 keV, and we 
used an interstellar H column density of 3.4 x 10^^ cm“^ 
which was obtai ned from Swift data analysis of the source 
(jLin et al.ll2013[ l. 

We obtained the background spectra for the bursts 
from pre-burst data and, when possible, also from post¬ 
burst data. Eor the tails we extracted the background 
spectra only from the pre-burst data. In all cases we 
used all layers of operating PGUs and hnally grouped 
the burst and tail spectra such that each spectral bin 
would contain at least 20 counts (except the burst in 
event B and burst 1 in event D, which have less counts 
than other bursts; these are grouped to contain at least 
10 counts). 

The durations of the four events with tails A, B, C, and 
D are 15, 3534, 624 and 773 s, respectively, determined 
based on Bayesian block representation. As mentioned in 
the previous section, the reason for the long tail duration 
of event B is that the algorithm could not find a block 
that goes below the assigned background level before the 
end of RXTE orbit. We note, however, that the emission 
after ~200 s is quite weak; in fact, inclusion of data after 
^200 s in the spectral analysis did not signihcantly alter 
the parameters. Therefore, we limit our investigation by 
that time. 

For each of the four events, we analyzed the burst and 
tail spectra separately. In the cases of events B, G, and D, 
the sufficiently long durations of the tails enabled us to 
investigate the spectral evolutions throughout the tails. 
To this end, we divided the tails of B and D into two 
segments, and tail of event G into three segments, sepa¬ 
rated by the other two bursts in this tail (see Figures |3l 
|4|and|5]), and performed time-resolved spectral analysis. 
We present the durations and count rates of the burst 
and tail spectra in Table [2] 

For modeling the burst and tail spectra, we employed 
a set of spectral models that are commonly used for mag- 
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Fig. 2.— Top panel: 2—20 keV light curves of the four events with 1 s resolution. The four events are identified here with labels A—D and 
indicated in each panel. Bottom panel: Bayesian block representation of the light curve. The blue dashed line represents the background 
level. Red vertical dotted lines indicate the start and end times of the events. For event B the tail end time is the end of the data (see 
text). 


netar spectral analysis: single blackbody, double black- 
body, optically thin thermal bremsstrahlung, power law, 
cutoff po wer law, as we l l as combinations of these model s 
(see, e.g.. Ilsrael et alJ (I2()10D: Ivon Kienlin et alJ (|2012D : 
Ivan der Horst et al.l (|2012D : ILin et al.l (|2012l) '). Below, we 
describe the spectral analysis results of each event. We 
also present the best-fit spectral parameters determined 
by the statistics in Table [3] and spectral evolution of 
events B, C, and D in Figures EHH and (5] We note that 
uncertainties are reported at the Icr level throughout the 


paper. 

4.1. Results of Spectral Analysis 
4.1.1. Event A 

The burst in the beginning of event A was saturated 
due to high number of burst photons. We therefore 
excluded the time intervals (a total of ^0.15 s) during 
which the count rate exceeded 18,000 c s“^PCU“^. The 
burst spectrum is described best with the cutoff power- 
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law model with a photon index of 0.68±0.08 and cutoff 
energy 14.24^11^ ixled = 1-14 with 49 degrees of 
freedom (dof)). We found that the tail of this event is 
fitted with a power-law model the best: index of 
(Xred = with 34 dof). The 2—25 keV unabsorbed 
fluxes of the burst and the tail are >9.87±0.12x 10“® 
and 9.23±0.44 x 10“^° erg cm“^ s“^, respectively. The 
corresponding energies, assumin g isotropic emission and 
the source distance of 5 kpc (iTiengo et al.l I2010H are 
>5.39x10^® erg and 4.69x10®^ erg for the burst and the 
tail, respectively. Note that the burst energy is a lower 
bound since the detector was saturated, and its true en¬ 
ergy content was higher than our estimate. 

4.1.2. Event B 

The burst in event B can be fitted with a single- 
blackbody model with a temperature of 4.86 iI1q 7 ® keV 
(Xred = 0-43 with 6 dof). The time-integrated spectrum 
of the tail of this event is also best fitted with a black- 
body model of kT = 2.22to'ig keV {xled — with 
53 dof). Note that the power-law model for this time- 
integrated tail yields a similar fit statistics: Xred 1-01 
with an index of 1.51±0.18. We also analyzed the spectra 
of each segment of the tail and found that both segments 
are described well with a blackbody model of changing 
temperature of d.Ol^Q^g and keV (see Figure 

ED- The corresponding radii of the blackbody emitting 
region are 0.25to;ii km and 0.36to;22 respectively 
(given d = 5 kpc). 

4.1.3. Event C 

Event C contains three bursts, one in the beginning of 
the event, the second one ~150 s after the first burst, and 
the last one separated by ^190 s from the second burst 
(see Figure S]). The first burst has also saturated the de¬ 
tector due to a high number of incoming burst photons, 
thus, we excluded the saturated portion of the burst (a 
total of ~0.18 s), applying the same count rate criterion 
used for the event A. For this burst, a combination of 
blackbody and power-law models provided the best fit 
(Xred = 1-24 with 52 dof); the blackbody temperature is 
15.46^2 58 power-law index is 2.01±0.14. This 

corresponds to a blackbody radius of 0 . 641 q; 4 q km. The 
slightly large Xred obtained in this best fit was actually 
due to the data around 14 keV. Excluding the energy 
range of 11—16 keV from the spectral analysis improves 
the fit statistics. We will investigate this burst in detail in 
§a§maz Mu§, S. et al. (2015, in preparation). Second and 
third bursts can be fitted well with a power-law model 
with photon indices 1.43±0.11 (x^ed = with 13 dof) 
and 1.79±0.16 (x^ed = with 8 dof), respectively. As 
for the tail, a blackbody model fits the time-integrated 
tail spectrum very well (Xred = 1-18 with 54 dof). The 
blackbody temperature is 2.12to;iJ keV, and correspond¬ 
ing radius is 0.24^® km. The three segments of this tail 
can also be fitted with the blackbody model, resulting in 
temperatures of 2.37±0.11, 2.17to;i5 and l.GOlQ ^g keV, 
respectively (see Figure |4]). These correspond to black¬ 
body emitting region radii of 0.291°' 44 , 0.241°!® and 
0.25ll!^ km. 


4.1.4. Event D 

At the onset of event D, there are two bursts sepa¬ 
rated by ^11 s. The first burst is best fitted with a 
single blackbody with a temperature of 5.13llg° keV 
but still resulted in large x^ (x^d = with 4 dof). 
Further investigation of this burst revealed a possible 
spectral feature around 14 keV. Similar spectral fea¬ 
tures around the same energy have been observed in 
other magnetars (iGavriil et al.ll200^ iWoods et al.l 120051 : 
iGavriil et al.ll2006l:IAn et al.ll2014ll . Inclusion of a Gaus¬ 
sian line improves the fit statistics significantly, but en¬ 
ergetics are comparable with or without the Gaussian 
line. This will also be investigated in detail in §a§maz 
Mu§, S. et al. (2015, in preparation). The second burst 
of this event is fitted with a single blackbody with kT 
= 7 . 53 ! 4 ; 3 g (Xred = 0-34 with 5 dof). The time- 
integrated tail spectrum is also fitted well with a black¬ 
body model with temperature 2 . 59 !g !2 keV (x?e, = 0.91 
with 53 dof). The corresponding blackbody radius is cal¬ 
culated as O.lStg!® km. Similar to the other events, 
we found that the tail can be modeled with a black¬ 
body of decreasing temperature, 3.23!g!3 keV in the 
first and 2 . 20 !°!g keV in the second segment (see Fig¬ 
ure 0 . The corresponding blackbody radii are 0.22±0.08 
and 0.19!g!g km, respectively. 

5. TEMPORAL ANALYSIS 
5.1. Phases of Bursts 

To investigate whether bursts that exhibit extended 
tails have any dependence on the spin phase of the neu¬ 
tron star, we have calculated the corresponding phases 
of the four particular bursts studied here. For this pur¬ 
pose, we used contempora.neou s phase-connected spin 
ephemerides of iDib et ahl l)2012f) for the events A and 
B, and spin ephem erides or frequencies provided by 
iKuiper et al.1 (I2012D for the remaining two events. We 
find that the energetic burst of event A starts at the 
spin phase, (j) of 0.25, which is the beginning of the 
peak plateau of the pulse profile (see Figure 0, continues 
throughout the pulse peak, and declines rapidly after f) of 
0.75, which is the end of the peak plateau. Event B spans 
the spin phase interval of 0.60 to 0.64. Event C starts at 
(j) = 0 . 11 , which is the rising part of the pulse profile, and 
ends somewhere during the peak of the persistent emis¬ 
sion. Finally, event D has two spikes separated by 11 s: 
the first one spans between (j) of 0.91 and 0.95, which is 
within the minimum phase of the pulse profile, and the 
second one spans between ef = 0.10 and 0 . 22 , which is 
again during the rising portion of the pulse profile. 

5.2. Dependence of Pulse Amplitude on Burst Ignition 

Energetic bursts preceding extended tails may induce 
observable changes in the pulse properties of the neutron 
star. In order to investigate this for SGR J1550—5418, we 
compared the pulse amplitudes prior to the onset of each 
of the four events and of those during the tails follow¬ 
ing the bursts. We created the phase-folded light curves 
in the 2—20 keV energy range considering only PCU 2 
for the entire tail and the pre-burst emission of each of 
the four bursts u sing the spin ephemeride s provided by 
iDib et ail l)2012ll and iKuioer et al.l (|2012f) accordingly. 


























6 


TABLE 3 

Best-fit Spectral Parameters (2—25 keV) for the Bursts and the Tails 


Event 

Model® 

kT 

(keV) 

Index 

BB Radius 
(km) 

Flux (erg cm ^ s 
(2—25 keV) Unabsorbed 

Isotropic Energy (erg) 

(2-10 keV) (2-25 keV) 


A 

Burst 

Cutoff PL*' 

- 

0.68+0.08 

- 

9.87+0.12 X 10-® 

2.43+0.03 X 10°® 

5.40+0.07 X 10°® 

1.14(49) 

Tail 

PL 

— 

1 q' 7 + 0.10 
-‘-•'^'-0.11 

— 

9.23+0.44 X 10-1° 

2.11+0.16 X 10°7 

4.69to 22 X 10°7 

0.84(34) 

B 

Burst 

BB 

4.86lJ;“« 

- 

-0.43 

8.391 j X 10"'’ 

7A2t\il X 10°® 

3.0ll°;i« X 10°° 

0.43(6) 

Tail 1-1-2 

BB 

r\ Qr)-|-0.22 

^•^^-0.19 

- 

0.28l“;== 

8.06l°;™ X 10-11 

3.17+0.31 X 10°7 

4.68l°;i° X 10°7 

0.98(53) 

Tail 1 

BB 

'^•^-‘--0.18 

- 

0.25l°+ 

2.14l°l® X 10-1° 

- 

- 

0.91(98) 

Tail 2 

BB 

-| yQ + 0.22 

-*-■ ''^-0.19 

— 

U.OD q 22 

4.7ll°;g X 10-11 

— 

— 

0.91(98) 

C 

Burst 1 

BB+PL 

15AGttll 

2.01+0.14 

U.04t_Q 4Q 

1.88l° °8 X 10-7 

6.12+0.16 X 10°7 

1.97t°;°7 X 10°® 

1.24(52) 

Burst 2 

PL 

- 

1.43+0.11 

- 

1.89+0.11 X 10-® 

4.50+0.33 X 10°° 

S-Olto.M X 10°° 

1.01(13) 

Burst 3 

PL 

- 

1.79+0.16 

- 

1.08+0.07 X 10-® 

3.70j:°;3° X 10°® 

6.45+0.43 X 10°® 

1.00(8) 

Tail 1+2+3 

BB 

9 1 9+0.11 
^•^^-0.10 

- 

n 

O.^4_0.io 

s.oelo jf X 10-11 

6.60+0.35 X 10°7 

9.32l°;°J X 10°7 

1.18(54) 

Tail 1 

BB 

2.37+0.11 

- 

o-^y_o.ii 

l-12lo.05 X 10-^° 

- 

- 

1.18(162) 

Tail 2 

BB 

9 1 '7+0.16 

^•^'-0.15 

- 

G.^^_0.11 

5.18l°;iJ X 10-11 

- 

- 

1.18(162) 

Tail 3 

BB 

i.Dy o ig 

— 

U.ZO Q 15 

2.07l°;29 X 10-11 

— 

— 

1.18(162) 

D 

Burst 1 

BB 

tr 1 q+1.50 
'J--L'J_o.94 

- 

0 43+0-37 
^•^'3_o.3o 

4.14 1°'®° X 10-® 

8.53+1.80 X 10°° 

3.7ltg;gg X 10°° 

1.56(4) 

Burst 2 

BB 

7 

< 3g 

- 

0.39l°;“ 

9.63l°|i X 10-° 

l-OOto.l® X 10°® 

1.30+0.13 X 10°7 

0.34(5) 

Tail 1+2 

BB 

z.oy_Q 

- 

o-isin? 

6.03l°;3i X 10-11 

7.83+0.42 X 10°7 

1.37+0.07 X 10°® 

0.91(53) 

Tail 1 

BB 

q 9q+0.14 
o.zo_Q 13 

- 

0.22+0.08 

2.22+0.11 X 10-1° 

- 

- 

1.12(103) 

Tail 2 

BB 

^■^^-0.18 

— 

01Q+0+^ 

u.iy o.io 

3.56l°;|° X 10-11 

— 

— 

1.12(103) 


® PL: power-law model; BB: blackbody model. 

Cutoff energy of this model is 14.24+j |® keV. 

We present the simultaneous fit results for the tail segments. 


We modeled the folded profiles using a sinusoidal model 
with a harmonic and computed the rms fractional am¬ 
plitude from the profiles. For event A, the fractional 
rms amplitude before the burst and during the tail were 
1.58±0.6% and 17.66±2.73%, respectively. For event B 
the pre-burst emission available in the same RXTE orbit 
was very short, and we estimated the pre-burst pulsa¬ 
tion using the data taken during the previous orbit. For 
this burst, we hnd that the pulsation amplitudes during 
pre-burst and tail were 2.32±1.15% and 13.68±2.38%, re¬ 
spectively. For event C, we obtained the pulse amplitude 
before the burst as 6.51±2.50% and in the tail it became 
5.07±1.39%. For the final event D, the amplitudes in 
the pre-burst emission and the tail were 8.16±2.25% and 
12.20±1.75%, respectively. 

It was shown for SGR 1900-1-14 that the pulse frac¬ 
tion increases immediately after t he burst and varie s 
with time in the tail of the burst (|Lenters et al.l[200^ . 
In order to check whether the rms amplitude of pulsa¬ 
tion varies during the tail itself for SGR J1550—5418, we 


considered smaller time segments during the tail of all 
the bursts except the first one and computed the phase- 
folded profile and rms fractional amplitudes in each of the 
short time segments. For event D, the rms amplitude is 
significantly high (~0.25) at the very beginning of the 
tail just following the burst and then quickly decreases 
(see Figure [7]). The high rms amplitude in the first 50 
s of the tail complies with the signif i cant detection of 
high-power pulsation bv lKuiner et al.l (|2012f) in the first 
52 s since the onset of that particular burst. After the 
first 100 s, the fractional rms amplitude value and its 
variation are consistent with those observed for the pre¬ 
burst emission. For event B the average rms amplitude 
showed a marginal increase during the early episodes of 
the tail compared to its pre-burst emission, though the 
increase was not as significant as in the case of event D. 
For event C, the rms fractional amplitude of pulsation 
did not vary significantly, remaining generally consistent 
with the variation of pulse amplitude during the pre¬ 
burst emission. To give a measure of the variation of 
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Fig. 3. — Spectral parameters of the burst and tail of event B. 
Top panel: 1 s binned light curve of the event. Dotted vertical lines 
indicate the time intervals of tail segments. Second panel: black- 
body temperature of the burst (red diamond) and tail segments 
(blue crosses). Third panel: blackbody radii of the burst and tail 
segments. Bottom panel: unabsorbed fluxes of the burst and tail 
segments. 

the amplitude of pulsation, we quote the mean and stan- 
dard deviation of the pulse amplitudes calculated over 
the finer time segments. The average fractional rms am¬ 
plitudes during pre-burst and the entire duration of tails 
for the three events are: 6.24±4.04% and 16.68±1.26% 
for event B, 7.93±4.03% and 6.90±3.39% for event C, 
and 9.64±2.80% and 13.25±6.32% for event D. Note that 
the errors are calculated from the deviations of ampli¬ 
tudes over the short time segments, over which the pulsed 
signal and fractional rms amplitude could only be poorly 
constrained. 

5.3. Search for QPOs in the Tails of Bursts 

We investigated the power spectra during the tails of 
the four events to search for quasi-periodic oscillations 
(QPOs) that have been first observed from SGR 1806—20 
and SGR 1900-1-14 among i nagnetars (llsrael et alJl200^ 
iStrohmaver fc Watts! 1200^ iWatts fc Strohmaved l2006ll . 
For this purpose, we computed the Leahy-normalized 
power spectra over each 1 s segment during the tail in the 
energy ranges 2—10 and 10—30 keV. For QPO detection, 
we chose a threshold of 3cr significance level, considering 
the number of trials as the number of frequency bins and 
the number of power spectra searched. Using our crite¬ 
rion, we did not detect any significant QPO features in 
the individual or averaged power spectra during the burst 


Fig. 4.— Spectral parameters of the bursts and tail of event C. 
Dotted vertical lines indicate the time intervals of tail segments. 
Top panel: 1 s binned light curve of the event. Second panel: 
blackbody temperature of the tail segments. Third panel: black- 
body radii of the tail segments. Fourth panel: unabsorbed fluxes of 
the bursts. Bottom panel: unabsorbed fluxes of the tail segments. 

tails, and the power spectra were consistent with Poisson 
noise. To further quantify the temporal behavior of the 
burst tails, we calculated the total noise rms power over 
the energy bands 2—10 and 10—30 keV. We divided the 
burst tails in 15 temporal segments and averaged over 
them to get the resultant power spectra in each band. 
We obtained the total rms noise in the 10—1020 Hz in¬ 
terval for the two energy bands for event A as 51.09% and 
51.05%, respectively; for event B as 53.28% and 53.11%, 
respectively; for event C as 54.24% and 53.24%, respec¬ 
tively; and finally for event D as 51.96% and 51.46%, 
respectively. We also computed the confidence level and 
corresponding chance occurrence probability of the high¬ 
est power detected in the power spectrum of each tail, 
again in the two energy bands. As we present in Ta¬ 
ble m there is no indication of periodic or quasi-periodic 
processes in the tail of these events. 

6. DISCUSSION 

We have searched for extended burst tails using the 
Bayesian blocks technique in a large collection of RXTE 
data of SGR J1550—5418. We identified four events with 
durations of 15, 3534, 624, and 773 s, and we have stud¬ 
ied spectral and temporal properties of these four events. 
We compare below general properties of these events with 
one another, as well as those of extended burst tails de¬ 
tected from other magnetars, namely, SGR 1900-1-14 and 




















































TABLE 4 

Significance levels of the maximum power observed in the power spectra in the 10—1020 Hz frequency range. 


Event 

Total rms 

Chance 

Number of 

Confidence 

Significance 


Noise Power 

Probability 

Trials 

Level (%) 

(o-) 


Energy Range 2—10 keV 

A 

51.09 

1.49x10"^ 

1144 

18.19 

0.23 

B 

53.28 

1.01x10-5 

13061 

87.65 

1.54 

C 

54.24 

9.44x10-5 

38924 

2.54 

0.03 

D 

51.96 

7.64x10-’^ 

51220 

96.16 

2.07 


Energy Range 10—30 keV 

A 

51.05 

3.58x10-^ 

1144 

66.38 

0.96 

B 

53.11 

5.74x10-5 

13061 

47.26 

0.63 

C 

53.24 

2.91x10-5 

38924 

32.21 

0.41 

D 

51.46 

1.98x10-5 

51220 

36.27 

0.47 



0 200 400 600 800 


Time (seconds since 505861869.378 s) 

Fig. 5.— Spectral parameters of the bursts and tail of event D. 
Dotted vertical lines indicate the time intervals of tail segments. 
Top panel: 1 s binned light curve of the event. Second panel: black- 
body temperature of the bursts (red diamonds) and tail segments 
(blue crosses). Third panel: blackbody radii of the bursts and tail 
segments. Fourth panel: unabsorbed fluxes of the bursts. Bottom 
panel: unabsorbed fluxes of tail segments. 

SGR 1806-20. 

The tail of event A has a much shorter duration 
(^13 s). Moreover, unlike events B, C and D, its X- 
ray spectrum has non-thermal character; it is described 
the best with a power-law model of index 1.37, resem¬ 
bling the SGR burst spectral shape in the RXTE/VGA 
passband. The other three burst tails have thermal 
spectra, described with a blackbody model, at temper¬ 
atures around 2—3 keV, and more importantly show¬ 
ing a clear trend of declining blackbody temperature 
over the course of the tail. According to the magne- 



Phase (cycles) 

Fig. 6.— Top panel: light curve and pulse profile of the tail 
detected in event A. Black bars in the pulse profile indicate the 
phase interval of the other events. Di and D 2 denote the two 
bursts observed in event D. 

tar model, the repeated bursts from SGR sources are 
likely due to fract uring of the solid neutron s tar crust 
by magnetic stress (iThomnson &: puncanlll995l) or mag¬ 
netic reconnection (|Lvutikovll2003fl . In either case, ener¬ 
getic electron—positron pair plasma would be induced 
into the magnetosphere which would be observed as 
energetic bursts. Note the important fact that event 
A was detected on 2009 January 22, that is the day 
SGR J1550—5418 was most burst active, and many en- 
erge tic and relatively long b ursts were also detecte d 
(see iMereghetti et al.l (I2009D . I Savchenko et al.l (I2010H ~). 
Therefore, event A is probably an exceptionally long 
burst, and the identified tail is most likely the contin¬ 
ued emission of this prolonged event. 

Extended burst tails were also seen from other mag- 
netars: two from SGR 1900-1-14, soon after its giant 
flare i n 1998 August and i ntermediate f l are in 2001 
April (|Ibrahim et al.l I200H iLenters et al.l 1200311 . and 
two from SGR 1806—20 during the 2003—2004 burst- 
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Time (seconds since 505861961.698741 s) 


Fig. 7.— Variation of rms fractional amplitude during the tail of 
event D. The rms fractional amplitudes were calculated at each 50 
s time bin during the tail. The dotted vertical line represents the 
separation between the tail segments considered for the spectral 
analysis. 

acti ve episode prior to its giant flare in 2004 Decem¬ 
ber (iGogiis et al.ll20lill . Energies of the bursts leading 
to extended tails in SGR 1900-1-14 were at the highest 
of the scale for typical magnetar bursts. On the other 
hand, preceding bursts of SGR 1806—20 tails were not 
the highest-energy ones: there were a lot of other ener¬ 
getic bursts that were not followed by tails. Similarly, we 
found that the energy of the tail triggering burst in SGR 
J1550—5418 could be as low as 3.01x1 0^^ erg, and there 
were a lot of more energetic bursts (jMereghetti et all 
120091 : Ivan der Horst et al.ll2012| l. The extended tail phe¬ 
nomenon is, therefore, a special case in magnetar bursts; 
there should be a minimum energy injection to ignite a 
tail, but not all bursts that are more energetic than the 
threshold would lead to extended tails. Besides, unlike 
the other magnetars, in the case of SGR J1550—5418 the 
ratios of total energy contained in the tails and in the 
bursts seem to vary from event to event. Hence, the ques¬ 
tion of what ignites extended tails cannot be answered 
simply with th e energetics of the m ain burst, as already 
pointed out by iGoeiis et al.l 1)2011!) . 
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An important common spectral property of all ex¬ 
tended tails observed now from three sources is that they 
all exhibit thermal spectra and the temperatures decline, 
that is, cooling throughout the tail. The spectral cool¬ 
ing behavior in SGR J1550—5418 is not as significant as 
in the cases of SGR 1900-1-14 and SGR 1806—20 tails, 
but still evident. The corresponding blackbody emit¬ 
ting area of SGR J1550—5418 tails remains fairly con¬ 
stant around 0.2—0.3 km for events B, G, and D. Similar 
to the earlier tails detected, SGR J1550—5418 extended 
tails are also exhibiting the cooling of a heated portion 
of the neutron star crust. Given the energetics argument 
above and the spectral nature of tails, we suggest that 
the source of heating is the bombardment of the neutron 
star surface with returning pairs in the trapped fireball 
(iThompson fc Duncanlll99^ . which could not efficiently 
radiate away. Therefore, the consequential cause of the 
generation of extended tails would depend on how effi¬ 
ciently the trapped plasma in the magnetosphere radiates 
away. 

As far as the temporal investigations of the extended 
tails are concerned, we do not find any phase depen¬ 
dence of leading bursts, and also any evidence of QPOs 
in the high-frequency domain, which are observed from 
the giant flares of SGR 1900-1-14 and SGR 1806—20, and 
are attrib uted to torsional oscillations of the neutron 
star crust (llsrael et al.ll2005l : iStrohmaver &: Wattd 120051 : 

IWatts &: Strohmaverl I2006lb Apart from other causes, 

the QPO phenomenon likely requires a much higher burst 
energetic threshold than that of extended tails to ex¬ 
cite torsional modes. Interestingly, we find a significant 
pulsed amplitude increase in the very early portion of 
event D, which has one of the least burst energetics. This 
burst energy argument conditionally strengthens our hy¬ 
pothesis of surface heating with returning pairs because 
as more energy is imparted back to the system, the en¬ 
ergy content of radiated portion would naturally be less. 
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